Determination of continuous variable entanglement by purity measurements 



Gerardo Adesso, Alessio Serafini, and Fabrizio Illuminati 

Dipartimento di Fisica "E. R. Caianiello", Universita di Salerno, INFM UdR di Salerno, 
INFN Sezione di Napoli, Gruppo Collegato di Salerno, Via S. Allende, 84081 Baronissi (SA), Italy 

(Dated: February 25, 2004) 

We classify the entanglement of two-mode Gaussian states according to their degree of total and partial 
mixedness. We derive exact bounds that determine maximally and minimally entangled states for fixed global 
and marginal purities. This characterization allows for an experimentally reliable estimate of continuous variable 
entanglement based on measurements of purity. 
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Quantum entanglement of Gaussian states constitutes a fun- 
damental resource in continuous variable (CV) quantum infor- 
mation [1]. Therefore, the quest for a theoretically satisfying 
and experimentally realizable quantification of the entangle- 
ment for such states stands as a major issue in the field. On the 
theoretical ground, a proper, computable quantitative charac- 
terization of the entanglement of Gaussian states is provided 
by the logarithmic negativity |2||. Experimental schemes to 
determine the entanglement of Gaussian states have been pro- 
posed both in the two-mode [3] and in the multipartite [4] 
instance. However, these schemes are based on homodyne 
detections, and require a full reconstruction of all the second 
moments of the Gaussian field. 

In this work, we present a theoretical framework to esti- 
mate the entanglement of two-mode Gaussian states by the 
knowledge of the total and of the two partial purities. This 
is achieved by deriving analytical a priori upper and lower 
bounds on the logarithmic negativity for fixed values of the 
global and marginal purities. We then show that the set of en- 
tangled Gaussian states is tightly contained between two ex- 
tremal surfaces of maximally and minimally entangled states. 
This quantification allows for a simple strategy to measure 
the entanglement of Gaussian states with reliable experimen- 
tal accuracy. In fact, measurements of global and marginal 
purities do not require the demanding reconstruction of the 
full covariance matrix and can be performed directly by ex- 
ploiting the technology of quantum networks 1 5 ] . 

Let us consider a two-mode continuous variable system, 
described by the Hilbert space Ti = Hi <£> TL2 resulting from 
the tensor product of the Fock spaces 7i k s. We denote by a k 
the annihilation operator acting on the space Tt k . Likewise, 
x k = (ofc + a\)/y/2 and p k = ~i(a k - %)/V2 are the 
quadrature phase operators of the mode k, the corresponding 
phase space variables being x k andp k . 

In the following, we will make use of the Wigner quasi- 
probability representation W(xi,Pi), defined as the Fourier 
transform of the symmetrically ordered characteristic func- 
tion. In Wigner phase space picture, the tensor product 7i = 
Ti-i (S> Ti-2 results in the direct sum V = T\ © T2 of the related 
phase spaces IYs. A symplectic transformation acting on the 
global phase space V corresponds to a unitary operator acting 
on H 1 6]. We will refer to a transformation Si = 5i © £2, 
with each Si G Spt2,R) acting on Ti, as to a "local symplectic 



operation", corresponding to a "local unitary transformation" 
Ui = U\®U2- The set of Gaussian states is defined as the set 
of states with Gaussian Wigner function 



W(X) = 



,-iJte 



'X 1 



7rvDet cr 



(1) 



where X = (x\ ,pi,X2,P2) € T, and we will denote by X the 
vector of operators (xi , p\, X2, pz). First moments have been 
neglected, since they can be set to zero by means of a local 
unitary transformation. Second moments form the covariance 
matrix cr of the Gaussian state er^ = h{XiXj + XjXi) — 
(Xi)(Xj). For simplicity, in what follows cr will refer both 
to the Gaussian state and to its covariance matrix. It is conve- 
nient to express cr in terms of the three 2x2 matrices a, (3, 
7 
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Heisenberg uncertainty principle can be expressed as [ 



-n > 
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where f2 = uj © u) is the usual symplectic form with Wjj = 
<5jj_i— <%+i, i, j = 1,2. Ineq. (fj) can be recast as a constraint 
on the SW/^r) invariants A = Dct a. + Dct (3 + 2 Det 7, and 
Deter 
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In general, the Wigner function transforms as a scalar under 
symplectic operations, while the covariance matrix cr trans- 
forms according to cr — »• S T crS, with S € Spurn. As it 
is well known |7], for any covariance matrix cr there exists a 
local canonical operation Si = Si © S2 that recasts cr in the 
"standard form" cr s f with a = diag {a, a}, (3 — diag {6, b}, 
7 = diag {c+,c_}, where a, b, c + , c_ are determined by 
the four local symplectic invariants Det cr, Det a, Det (3, and 
Det 7. 

Any bipartite Gaussian state cr can always be writ- 
ten as cr = S T uS for some S S Sp&jt and v = 
diag n_, n + , n + }. The quantities n T constitute the 



2 



symplectic spectrum of er; they are determined by the global 
symplectic invariants | 



2n% = A =F V A 2 -4 Deter . 



(5) 



In terms of Ineq. @ becomes simply n_ > 1/2. 

We will characterize the mixedness of a quantum state g 
by its purity [i = Trg 2 . For a rt-mode Gaussian state er 
the purity is simply evaluated integrating the Wigner function, 
yielding ^ = 1/(2" VDet a). 

As for the entanglement, we recall that the positivity of the 
partially transposed (PPT) state <x is equivalent to separability 
for any two-mode Gaussian state er [6]. In terms of symplec- 
tic invariants, partial transposition corresponds to flipping the 
sign of Det 7, so that A turns into A = A — 4 Det 7. The 
symplectic spectrum n T of cr is simply found inserting A for 
A in Eq. (0. If n_ is the smallest symplectic eigenvalue of 
the partially transposed covariance matrix e>, a state er is sep- 
arable if and only if 



> 1/2 . 



(6) 



A bona fide measure of entanglement for two-mode Gaus- 
sian states should thus be a monotonically decreasing function 
of n_, quantifying the violation of inequality 0. A com- 
putable entanglement monotone for generic two-mode Gaus- 
sian states is provided by the logarithmic negativity Ej\f = 
max{0, — ln(2n_)} [2]. For symmetric Gaussian states, i.e. 
states whose standard form is characterized by a = (3, an- 
other computable entanglement monotone is provided by the 
entanglement of formation 1 10]. However, in this subcase the 
two measures provide the same characterization of entangle- 
ment and are fully equivalent. Therefore, from now on we will 
adopt the logarithmic negativity to quantify the entanglement 
of two-mode Gaussian states. 

We now show that a generic state in standard form can be 
reparametrized in terms of the Sp^.R) invariants /1 (the global 
purity) and A, and of the Sp^M) © 5p(2,R) invariants ^1 and 
/12, where /i^ denotes the purity of the reduced state in mode 
i {i = 1,2). For a generic two-mode Gaussian state cr s f we 
thus have 
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Deter = (ab) 2 - ab(c 2 + +c 2 _) + (c+c_) 2 , (8) 

a 2 + b 2 + 2c+c_ . (9) 

Eqs. J7I9I are easily inverted to provide the following 
parametrization 
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with eS ' fc, 2 ) 2 -4^lA] 
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The global and marginal purities range from to 1, con- 
strained by the condition 



(12) 



a direct consequence of Heisenberg uncertainty relations. It 
implies that no Gaussian LPTP (less pure than product) states 
exist, at variance with the case of two-qubit systems fllll . 
Eqs. J5I10I ITTt determine the smallest symplectic eigenvalue 
of the covariance matrix er and of its partial transpose er 
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where A = — A + l/2/i 2 + l/2/i|. This parametrization 
describes physical states if the radicals in Eqs. J 1 1 1 1 1 3> exist 
and Ineq. ® , expressing the Heisenberg principle, is satisfied. 
All these conditions can be combined and recast as upper and 
lower bounds on the invariant A 

J_ (Mi - M2) 2 
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The invariant A has a direct physical interpretation: at given 
global and marginal purities, it determines the amount of en- 
tanglement of the state. In fact, one has 
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The smallest symplectic eigenvalue of the partially transposed 
state is strictly monotone in A. Therefore the entanglement of 
a generic Gaussian state er with global purity [i and marginal 
purities Mi,2 strictly increases with decreasing A. Since A 
has both lower and upper bounds, due to Ineq. J 1 4b . not only 
maximally but also minimally entangled Gaussian states exist. 
This is an important result concerning the relation between en- 
tanglement and purity of quantum states: the entanglement of 
a Gaussian state is tightly bound by the amount of global and 
marginal purities, with only one remaining degree of freedom 
related to the invariant A. 

We now aim to characterize extremal (maximally or mini- 
mally) entangled Gaussian states for fixed global and marginal 
purities. Let us first consider the states saturating the lower 
bound in Eq. (I14> . which entails maximal entanglement. They 
are Gaussian maximally entangled mixed states (GMEMS), 
admitting the following parametrization 
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We now recall that Gaussian squeezed thermal states are states 
of the form er = S r vS r , where S r is the symplectic repre- 
sentation of the two-mode squeezing operator Exp [r(a\a2 — 



3 



a\a\)/2], while v = diag{n_ 



.}. These states 



are in standard form with a = n_ cosh r + n + sinh r, b — 
n + cosh 2 r+n_ sinh 2 r, c + = — c_ = sinh2r/2. 
In the pure case {n T = 1/2) they reduce to two-mode 
squeezed vacua. We thus find that states of the form of 
Eq. d!6i are non-symmetrical squeezed thermal states with 
tanh2r = 0*iM2 — M1M2/M) Ami + Ma)- These states are 
separable in the range 



M< 



MlM2 



Ml + M2 - M1M2 
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In such a separable region in the space of purities, no en- 
tanglement can occur for states of the form of Eq. d!6l >. 
while, outside this region, they are GMEMS. We now con- 
sider the class of states that saturate the upper bound in 
Eq. HAl . They determine the class of Gaussian least entan- 
gled mixed states (GLEMS). Violation of Ineq. d!7l > implies 
that (1 + I//1 2 ) /4 < + /i 2 ) 2 / 4 MiM2 - W Therefore, 
outside the separable region, GLEMS fulfill 
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Eq. dl 81 expresses saturation of Heisenberg relation (|4j. We 
thus find that the most semiclassical states of minimum quan- 
tum uncertainty are Gaussian least entangled states. GLEMS 
in standard form are characterized by 
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According to the PPT criterion, GLEMS are separable only 
for,u < p 1 p 2 /\/ r i4' 
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both separable and entangled states can be found. The very 
narrow region defined by Ineq. d20i is the only coexistence re- 
gion between entangled and separable Gaussian mixed states. 
Furthermore, Ineq. d!4l> leads to the following constraint on 
the purities 
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For purities which saturate Ineq. lEIt . GMEMS and GLEMS 
coincide and we have a unique class of states depending 
only on the marginal purities Ml, 2- They are Gaussian max- 
imally entangled states for fixed marginals (GMEMMS). The 
maximal entanglement of a Gaussian state decreases rapidly 
with increasing difference of marginal purities, in analogy 
with finite-dimensional systems II ill . For symmetric states 
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FIG. 1: Summary of entanglement properties of symmetric Gaus- 
sian states with given global and marginal purities. In the entangled 
region, the average logarithmic negativity Eq. 1251 is depicted, grow- 
ing from gray to black. The dashed area is the coexistence region of 
separable and entangled states. 



0*i = M2) Ineq. (12 1 1 reduces to the trivial bound fi < 1 and 
GMEMMS reduce to pure two-mode squeezed states. We can 
summarize the previous results in the following scheme, clas- 
sifying all the two-mode Gaussian physical states according 
to their degree of global and marginal purities 
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Knowledge of the global and marginal purities thus accurately 
characterizes the entanglement of Gaussian states, providing 
strong sufficient conditions and analytical bounds. As we will 
now show, marginal and global purities allow also an accurate 
quantification of entanglement. Outside the separable region, 
GMEMS attain maximum logarithmic negativity E^f max 
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while, in the entangled region (see Eq. \22Y ). GLEMS acquire 
minimum logarithmic negativity E_\f min 
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FIG. 2: Upper and lower bounds on the logarithmic negativity as 
functions of the global and marginal purities of symmetric Gaussian 
states. The black (gray) surface represents GMEMS (GLEMS). 



Knowledge of A (i.e. of the full covariance matrix) would 
allow for an exact quantification of the entanglement. How- 
ever, we will now show that an estimate based only on 
the knowledge of the experimentally measurable global and 
marginal purities turns out to be quite accurate. We can in 
fact quantify the entanglement of Gaussian states with given 
global and marginal purities by the "average logarithmic neg- 
ativity" £V defined as 



max 



(25) 



We can then also define the relative error 8Em on Em as 

En max: 



<S£V0l,2,M) 



En (Ml,2) + ENmin(Hl,2, ^) 



(26) 



It is easily seen that this error decreases both with increasing 
global purity and decreasing marginal purities, i.e. with in- 
creasing entanglement. For ease of graphical display, let us 
consider the important case of symmetric Gaussian states, for 
which the reduction jii — ji2 = Hi occurs. 
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FIG. 3: The relative error 5 Em Eq. 1261 on the average logarithmic 
negativity as a function of the ratio [i/fii, plotted at /i = 0.5. 



Fig-fflshows the classification of the entanglement of sym- 
metric states depending on their global and marginal purities. 
Notice in particular the very narrow region of coexistence of 
separable and entangled states. In Fig.|2j EMmini^ii A 4 ) °f 
Eq. J24t and .E^/max (/•**> A 1 ) °f Eq. d23i are plotted versus fa 
and \i. In the case fi = 1 the upper and lower bounds correctly 
coincide, since for pure states the entanglement is completely 
quantified by the marginal purity. For mixed states this is not 
the case, but, as the plot shows, knowledge of the global and 
marginal purities strictly bounds the entanglement both from 
above and from below. 



The relative error SEj^if^i, m) given by Eq. i26\ is plot- 
ted in Fig. [5] as a function of the ratio fj,//j,i. It decays expo- 
nentially, and falls below 5% for /i > /ij. Thus detection of 
genuinely entangled states is always assured by this method, 
except at most for a small set of states with very weak entan- 
glement (states with Em < 1). Moreover, the accuracy is even 
greater in the general non-symmetric case /ii ^ /12, because 
the maximal entanglement decreases in such an instance. The 
above analysis demonstrates that the average logarithmic neg- 
ativity Em is a reliable estimate of the logarithmic negativity 
Em, improving as the entanglement increases. This allows 
for an accurate quantification of CV entanglement by knowl- 
edge of the global and marginal purities. The latter quantities 
may be in turn amenable to direct experimental determination 
by exploiting the technology of quantum networks [5], even 
without homodyningl 12]. The present work thus may provide 
a powerful operative characterization and quantification of the 
entanglement of generic Gaussian states. 
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